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Abstract: We deal with the test of the general relativistic gravitomagnetic Lense-Thirring effect currently being 

conducted in the Earth's gravitational field with the combined nodes Q of the laser-ranged geodetic satellites 
LAGEOS and LAGEOS II. One of the most important sources of systematic uncertainty on the orbits of 
the LAGEOS satellites, with respect to the Lense-Thirring signature, is the bias due to the even zonal 
harmonic coefficients Ji of the multipolar expansion of the Earth's geopotential which account for the 
departures from sphericity of the terrestrial gravitational potential induced by the centrifugal effects of its 
diurnal rotation. The issue addressed here is: are the so far published evaluations of such a systematic 
error reliable and realistic? The answer is negative. Indeed, if the difference AJi among the even zonals 
estimated in different global solutions (EIGEN-GRACE02S, EIGEN-CG03C, GGM02S, GGM03S, ITG- 
Grace02, ITG-Grace03s, JEM01-RL03B, EGM2008, AIUB-GRACEOIS) is assumed for the uncertainties 
SJi instead of using their more-or-less calibrated covariances crj^, it turns out that the systematic error 5fi 
in the Lense-Thirring measurement is about 3 to 4 times larger than in the evaluations so far published 
based on the use of the covariances of one model at a time separately, amounting up to 37% for the pair 
EIGEN-GRACE02S/ITG-Grace03s. The comparison among the other recent GRACE-based models yields 
bias as large as about 25 — 30%. The major discrepancies still occur for J4, Jg and Jg, which are just to 
which the zonals the combined LAGEOS/LAGOES II nodes are most sensitive. 
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1. Introduction 

In the weak-field and slow motion approximation, tlie Einstein field equations of general relativity get linearized 
to a form resembling Maxwell's equations of electromagnetism. Thus, a gravitomagnetic field, induced by the 
off-diagonal components goi,i ~ 1,2,3 of the space-time metric tensor related to the mass-energy currents of 
the source of the gravitational field, arises [1]. It affects in several ways the motion of, e.g., test particles and 
electromagnetic waves [2]. Perhaps the most famous gravitomagnetic effects are gyroscope precession [3, 4] and 
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the Lense-Thirring^ precessions [6] of the orbit of a test particle, both occurring in tlie field of a central slowly 
rotating mass like a planet. 

The measurement of gyroscope precession in the Earth's gravitational field has been the goal of the dedicated 
space-based GP-B mission'^ [7, 8] launched in 2004; its data analysis is still in progress. 

In this paper we critically discuss some issues concerning the test of the Lense-Thirring effect performed with 
the LAGEOS and LAGEOS II terrestrial artificial satellites [9] tracked with the Satellite Laser Ranging (SLR) 
technique [10]. 

[11, 12] proposed measuring the Lense-Thirring nodal precession of a pair of counter-orbiting spacecraft in terres- 
trial polar orbits and equipped with drag-free apparatus. A somewhat equivalent, cheaper version of such an idea 
was put forth in Ref. [13] whose author proposed to launch a passive, geodetic satellite in an orbit identical to 
that of the LAGEOS satellite apart from the orbital planes which should have been displaced by 180 deg apart.'' 
The measurable quantity was, in this case, the sum of the nodes of LAGEOS and of the new spacecraft, later 
named LAGEOS III, LARES, WEBER-SAT, in order to cancel the confounding effects of the multipoles of the 
Newtonian part of the terrestrial gravitational potential (see below). Although extensively studied by various 
groups [14-16], such an idea has not been implemented for a long time. Recently, the Italian Space Agency (ASI) 
has approved this project and should launch a VEGA rocket for this purpose at the end of 2009-beginning of 2010 
(http: //www. asi . it /en/activity /cosmology /lares). For recent updates of the LARES/ WEBER-SAT mission, 
including recently added additional goals in fundamental physics and related criticisms, see Refs. [17-24]. 
Among scenarios involving existing orbiting bodies, the idea of measuring the Lense-Thirring node rate with the 
just launched LAGEOS satellite, along with the other SLR targets orbiting at that time, was proposed in Ref. 
[25]. Tests have been effectively performed using the LAGEOS and LAGEOS II satellites [26], according to a 
strategy [27] involving a suitable combination of the nodes of both satellites and the perigee lj of LAGEOS II. 
This was done to reduce the impact of the most relevant source of systematic bias, viz. the mismodelling in 
the even (I — 2,4, 6. . .) zonal (m — 0) harmonics Jf of the multipolar expansion of the Newtonian part of the 
terrestrial gravitational potential:* they account for non-sphericity of the terrestrial gravitational field induced by 
centrifugal effects of the Earth's diurnal rotation. The even zonals affect the node and the perigee of a terrestrial 
satellite with secular precessions which may mimic the Lense-Thirring signature. The three-elements combination 
used allowed for removing the uncertainties in J2 and J4. In [28] a ~ 20% test was reported by using the^ EGM96 

^ According to an interesting historical analysis recently performed in Ref. [5], it would be more correct to speak 
about an Einstein-Thirring-Lense effect. 

See http: // einstein. Stanford, edu/ 
^ LAGEOS was put into orbit in 1976, followed by its twin LAGEOS II in 1992. 

* The relation among the even zonals Ji and the normalized gravity coefficients Ceo is Ji — + 1 Cfo- 

^ Contrary to the subsequent CHAMP/GRACE-based models, EGM96 relies upon multidecadal tracking of SLR 
data of a constellation of geodetic satellites including LAGEOS and LAGEOS II as well; thus the possibility of a 
sort of a — priori 'imprinting' of the Lense-Thirring effect itself, not solved- for in EGM96, cannot be neglected. 
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[29] Earth gravity model; subsequent analyses showed that such an evaluation of the total error budget was 
overly optimistic in view of the likely unreliable computation of the total bias due to the even zonals [30-32] . An 
analogous, huge underestimation turned out to hold also for the effect of non-gravitational perturbations [33] like 
direct solar radiation pressure, the Earth's albedo, various subtle thermal effects depending on the the physical 
properties of the satellites' surfaces and their rotational state [31, 34-40], which the perigees of LAGEOS-like 
satellites are particularly sensitive to. As a result, the realistic total error budget in the test reported in Ref. [28] 
might be as large as 60 — 90% or even more (by considering EGM96 only). 

The observable used in Ref. [9] with the GRACE-only EIGEN-GRACE02S model [41] and in Ref. [42] with other 
global terrestrial gravity solutions was the following linear combination*' of the nodes of LAGEOS and LAGEOS 
II, explicitly computed in Ref. [44] following the approach proposed in Ref. [27]: 

r ALAGEOS , ALAGEOS II 

/ = S2 + CiS2 , (1) 



where 
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The coefficients of the aliasing classical node precessions [45] fidass = "Y^i^.iJi induced by even zonals have 
been analytically worked out in e.g. [30]; a,e,i are the satellite's semimajor axis, eccentricity and inclination, 
respectively and yield ci = 0.544 for eq. (2). The Lense-Thirring signature of eq. (1) amounts to 47.8 milliarc- 
seconds per year (mas yr^^). The combination eq. (1) allows, by construction, to remove the aliasing effects due 
to the static and time- varying parts of the first even zonal J2. The nominal bias (computed with the estimated 
values of Ji, £ = 4, 6...) due to the remaining higher degree even zonals would amount to about 10^ mas yr~^; the 
need of a careful and reliable modeling of such an important source of systematic bias is, thus, quite apparent. 
Conversely, the nodes of the LAGEOS-type spacecraft are affected by the non-gravitational accelerations ~ 1% 
of the Lense-Thirring effect [36-40] . For a comprehensive, up-to-date overview of the numerous and subtle issues 
concerning the measurement of the Lense-Thirring effect see [46]. 
Here, we will address the following questions: 

• Has the systematic error due to the competing secular node precessions induced by the static part of the 
even zonal harmonics been realistically evaluated so far in literature? (Section 2) 

• Are other approaches to extract the gravitomagnetic signature from the data feasible? (Section 3) 



See also [31, 32, 43]. 
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2. The systematic error of gravitational origin 

The realistic evaluation of the total error budget of such a test raised a lively debate [47-53], mainly focussed 
on the impact of the static and time-varying parts of the Newtonian component of the Earth's gravitational 
potential through the aliasing secular precessions induced on a satellite's node. A common feature of all the 
competing evaluations so far published is that the systematic bias due to the static component of the geopotential 
was calculated always by using the released (more or less accurately calibrated) covariances ctj^, of one Earth 
gravity model solution at a time for the uncertainties 5Ji in the even zonal harmonics, yielding a percentage error 
particular to each model. 

Since it is always difficult to reliably calibrate the formal, statistical uncertainties in the estimated zonals of the 
covariance matrix for a global solution, it is much more realistic and conservative to instead take the differences^ 
AJf between the estimated even zonals for different pairs of Earth gravity field solutions as representative of 
the real uncertainty SJi in the zonals [55]. In Table 1-Table 12 we present our results for the most recent 
GRACE-based models released so far by different institutions and retrievable on the Internet at* http:/ /legem. gfz- 
potsdam.de/ICGEM/ICGEM.html. The models used are EIGEN-GRACE02S [41] and EIGEN-CG03C [57] from 
GFZ (Potsdam, Germany), GGM02S [58] and GGM03S [56] from CSR (Austin, Texas), ITG-Grace02s [59] and 
ITG-Grace03 [60] from ICG (Bonn, Germany), JEM01-RL03B from JPL (NASA, USA), EGM2008 [61] from 
NGA (USA) and AIUB-GRACEOIS [62] from AIUB (Bern, Switzerland). This approach was taken also in Ref. 
[27] with the JGM3 and GEMT-2 models. We included both the sum of the absolute values (SAY) of each 
mismodelled term and the square root of the sum of the squares (RSS) of each mismodelled term. 
The systematic bias evaluated with a more realistic approach is about 3 to 4 times larger than one can obtain 
by only using this or that particular model. The scatter is still quite large and differs greatly from that 5 — 10% 
claimed in Ref. [9]. In particular, it appears that J4, Jg, and to a lesser extent Jg, the most relevant zonals for us 
owing to their effecta on the combination of eq. (1), are the most uncertain ones, with discrepancies AJ^ between 
different models generally larger than the sum of their covariances ctj^ whether calibrated or not. 
Our approach is valid also for all of the tests performed so far with the LAGEOS and LAGEOS II satellites. 
Another possible strategy, that takes into account the scatter among the various solutions, is to compute mean 
and standard deviation of the entire set of values of the even zonals for the models considered so far, degree by 
degree, and then to take the standard deviations as representative of the uncertainties 5,h,l = 4,6,8, .... This 
yields 5^ — 15%, slightly larger than that recently obtained in Ref. [42]. But in evaluating mean and standard 
deviation for each even zonals, the authors of Ref. [42] also used global gravity solutions like EIGEN-GL04C and 
EIGEN-GL05C which include data from the LAGEOS satellite itself; this may likely have introduced a sort of 

^ See Fig. 5 of [54] for a comparison of the estimated C40 m different models. 

* / thank J Ries, CSR, and M Watkins (JPL) for having provided me with the even zonals of the GGM03S [56] 
and JEM01-RL03B models. 
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Table 1. Impact of the mismodclling in the even zonal harmonics on ft ^ jo^^*^^*^^ + ciO^/^^^^ "| AJt, £ ^ 4, . . . , 20, in 
mas yr^^. Recall that Jn — —y^2i + 1 C^o; for the uncertainty in the even zonals we have taken here the difference 
ACio ^ - C^o ^ I between the model X ^ EIGEN-CG03C [57] and the model Y ^ EIGEN-GRACE02S [41]. EIGEN- 

CG03C combines data from CHAMP (860 days out of October 2000 to June 2003), GRACE (376 days out of February 
to May 2003, July to December 2003 and February to July 2004) and terrestrial measurements; EIGEN-GRACE02S is 
based on 110 days (out of August and November 2002 and April, May and August 2003) of GRACE-only GPS-GRACE 
high-low satellite-to-satcllite data, on-board measurements of non-gravitational accelerations, and especially GRACE 
intersatellite tracking data. CTx/y ^re the covariance calibrated errors for both models. Values of fi smaller than 0.1 
mas yr~^ have not been quoted. The Lense-Thirring precession of the combination of eq. (1) amounts to 47.8 mas 
yr~^. The percent bias Sfi has been computed by normalizing the linear sum of f^, £ — 4, . . . , 20 to the Lense-Thirring 
precession. The discrepancies between the models are significant since AC^o arc larger than the linearly added sigmas 
for £ ^ 4, ...16. 
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Table 2. Impact of the mismodclling in the even zonal harmonics as solved for in X— GGM02S [58] and Y— ITG-Grace02s [59]. 

GGM02S is based on 363 days of GRACE-only data (GPS and intersatellite tracking, neither constraints nor regular- 
ization applied) spread between April 4, 2002 and Dec 31, 2003. The a are formal for both models. Ado ai'c always 
larger than the linearly added sigmas, apart from i — 12 and i — 18. 
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favorable a priori "imprint" of the Lense-Thirring effect itself. Moreover, the authors of Ref. [42] gave only a 
RSS evaluation of the total bias. 

We must also remember to add the further bias due to the cross-coupling between J2 and the orbit inclination, 
evaluated to be about 9% in Ref. [52]. 
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Table 3. Impact of the mismodelhng in the even zonal harmonics as solved for in X— GGM02S [58] and Y— EIGEN-CG03C [57]. 

The (T are formal for GGM02S, calibrated for EIGEN-CG03C. AC^o arc always larger than the linearly added sigmas. 
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Table 4. Bias due to the mismodcUing in the even zonals of the models X— ITG-Grace03s [60], based on GRACE-only accumulated 
normal equations from data out of September 2002-April 2007 (neither apriori information nor regularization used), and 
Y— GGMQ2S [58]. The a for both models arc formal. AC^o are always larger than the linearly added sigmas, apart from 
l=^Vl and I ^ 18. 
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3. A new approach to extract the Lense-Thirring signature from 
the data 

The technique adopted so far by the authors of Ref. [9] and Ref. [42] to extract the gravitomagentic signal 
from the LAGEOS and LAGEOS II data is described in detail in Refs. [54, 63]. The Lense-Thirring force is 
not included in the dynamical force models used to fit the satellites' data. In the data reduction process no 
dedicated gravitomagnetic parameter is estimated, contrary to e.g. station coordinates, state vector, satellites' 
drag coefficients Cd and Cr, etc.; its effect is retrieved with a sort of post-post-fit analysis in which the time 
series of the computed^ "residuals" of the nodes with the difference between the orbital elements of consecutive 
arcs, combined with eq. (f), is fitted with a straight line. 

In order to enforce the reliability of the ongoing test it would be desirable to follow other approaches as well. For 



^ The expression "residuals of the nodes" is used, strictly speaking, in an improper sense because the Keplerian 
orbital elements are not directly measured quantities. 
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Table 5. Bias due to the mismodelling in the even zonals of the models X — GGM02S [58] and Y — GGM03S [56] retrieved from 
data spanning January 2003 to December 2006. The a for GGM03S are calibrated. AC^o a^re larger than the linearly 
added sigmas for £ — 4, 6. {The other zonals arc of no concern) 
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Table 6. Bias due to the mismodelling in the even zonals of the models X ^ EIGEN-GRACE02S [41] and Y ^ GGM03S [56]. The 
a for both models arc calibrated. AC^o ^^'^ always larger than the linearly added sigmas apart from I — 14, 18. 
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instance, the gravitomagnetic force could be modelled in terms of a dedicated solve-for parameter (not necessarily 
the usual PPN 7 one) which could be estimated in the least-squares sense along with all the other parameters 
usually determined, and the resulting correlations among them could be inspected. Or, one could consider the 
changes in the values of the complete set of the estimated parameters with and without the Lense-Thirring effect. 
A first, tentative step towards the implementation of a similar strategy with the LAGEOS satellites in term of 
the PPN parameter 7 has been recently taken in Ref. [64]. 



4. Conclusions 



In this paper we have shown how the so far published evaluations of the total systematic error in the Lense- 
Thirring measurement with the combined nodes of the LAGEOS satellites due to the classical node precessions 
induced by the even zonal harmonics of the geopotential are likely optimistic. Indeed, they are all based on the use 
of elements from the covariance matrix, more or less reliably calibrated, of various Earth gravity model solutions 
used one at a time separately in such a way that the model X yields an error of x%, the model Y yields an error 
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Table 7. Bias due to the mismodelling in the even zonals of the models X — JEM01-RL03B, based on 49 months of GRACE-only 
data, and Y — GGM03S [56]. The a for GGM03S arc calibrated. AC^o a^rc always larger than the linearly added sigmas 
apart from £ — 16. 
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Table 8. Bias due to the mismodelling in the even zonals of the models X ^ JEM01-RL03B and Y ^ ITG-Grace03s [60]. The cr 
for ITG-Grace03s arc formal. AC^o arc always larger than the linearly added sigmas. 
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y%, etc. Instead, comparing the estimated values of the even zonals for different pairs of models allows for a 
much more realistic evaluation of the real uncertainties in our knowledge of the static part of the geopotential. 
As a consequence, the bias in the Lense-Thirring effect measurement is about three to four times larger than that 
so far claimed, amounting to tens of parts per cent (37% for the pair EIGEN-GRACE02S and ITG-GRACE03s, 
about 25-30% for the other most recent GRACE- based solutions). 

Finally, we have pointed out the need of following different strategies in extracting the Lense-Thirring pattern 
from the data; for instance by explicitly modelling it in fitting the SLR data of LAGEOS and LAGEOS II, and 
estimating the associated solve-for parameter in a least-square sense along with the other parameters usually 
determined. 



References 



Lorenzo lorio 



Table 9. Aliasing effect of the mismodelling in the even zonal harmonics estimated in the X— ITG-Grace03s [60] and the Y— EIGEN- 
GRACE02S [41] models. The covariance matrix a for ITG-Grace03s arc formal, while the ones of EIGEN-GRACE02S 
are calibrated. AC^o ^^rc larger than the linearly added sigmas for ^ — 4, 20, apart from I — 18. 
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Table 10. impact of the mismodelling in the even zonal harmonics estimated in the X— EGM2008 [61] and the Y— EIGEN- 
GRACE02S [41] models. The covariance matrix a are calibrated for both EGM2008 and EIGEN-GRACE02S. AC^o 
are larger than the linearly added sigmas for £ — 4, 20, apart from I — 18. 



i ACeo 


(EGM2008-EIGEN-GRACE02S) 


0"X + cry 


fe (mas yr-l) 


4 


2.71 X 10-11 


8.3 X 10-12 


10.0 


6 


2.35 X 10-11 


4.1 X 10-12 


5.0 


8 


1.23 X 10-11 


2.7 X 10-12 


0.4 


10 


9.2 X 10-12 


2.9 X 10-12 


0.1 


12 


4.1 X 10-12 


1.9 X 10-12 




14 


5.8 X 10-12 


1.8 X 10-12 




16 


3.4 X 10-12 


1.5 X 10-12 




18 


5 X 10-13 


1.5 X 10-12 




20 


1.8 X 10-12 


1.5 X 10-12 






S^j. = 33% (SAV) 


Sfi = 23% (RSS) 
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Table 11. Bias due to the mismodclling in the even zonals of the models X — JEM01-RL03B, based on 49 months of GRACE-only 
data, and Y — AIUB-GR ACEOlS [62], The latter one was obtained from GPS satcllitc-to-satcllitc tracking data and 
K-band range-rate data out of the period January 2003 to December 2003 using the Celestial Mechanics Approach. 
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formal. AC^o a^rc always larger than the linearly added sigmas. 





{JEM01-RL03B-AIUB-GRACE01S) 


0"X + 0"Y 


ft (mas yr-i) 


4 


2.95 X IQ-ii 


2.1 X 10-12 


11 


6 


3.5 X 10-12 


1.3 X 10-12 


0.8 


8 


2.14 X 10"" 


5 X 10-13 


0.7 


10 


4.8 X 10-12 


5 X 10-13 




12 
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5 X 10-13 




14 
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5 X 10-13 




16 


8 X 10-13 
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18 


7 X 10-13 


5 X 10-13 




20 


1.0 X 10-12 


5 X 10-13 






SfjL = 26% (SAV) 


5fi = 23% (RSS) 





Table 12. Bias due to the mismodelling in the even zonals of the models X = EIGEN-GRACE02S [41] and Y = J^IUB-GRACEOIS 
[62]. The a for AIUB-GRACEOIS are formal, while those of EIGEN-GRACE02S are calibrated. MJ to are larger than 
the linearly added sigmas for £ — 4, 6, 8, 16. 
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2 X 10-13 


1.5 X 10-12 





5n = 34% (SAV) 5fi = 25% (RSS) 
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